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Time dependent friction in a free gas 

Cristiano Fanelli* Francesco Sisti, + Gabriele V. Stagno -t 


Abstract 

We consider a body immersed in a perfect gas, moving under the action of a constant force E 
along the x axis . We assume the gas to be described by the mean-field approximation and interacting 
elastically with the body, we study the friction exerted by the gas on the body fixed at constant velocities. 
The dynamic in this setting urns studied in £3}, @ and £5} for object with simple shape, the first study 
where a simple kind of concavity was considered was in £Z0}, showing new features in the dynamic but 
not in the friction term. The case of more general shape of the body was left out for further difficulties, 
we believe indeed that there are actually non trivial issues to be faced for these more general cases. 

To show this and in the spirit of getting a more realistic perspective in the study of friction problems, 
in this paper we focused our attention on the friction term itself, studying its behavior on a body with a 
more general kind of concavity and fixed at constant velocities. We derive the expression of the friction 
term for constant velocities, we show how it is time dependent and we give its exact estimate in time. Fi¬ 
nally we use this result to show the absence of a stationary velocity in the actual dynamic of such a body. 
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1. Introduction 

We consider a body moving at a constant velocity in in a homogeneous fluid at rest, under the 
action of a constant force E, The gas is initially at thermal equilibrium and it interacts through 
elastic collisions with the body. The aim is to show the effect of the body's shape to the friction 
acting on it, giving an exact estimate of it. 

In 0 a model of free gas of light particles elastically interacting with a simple shaped 
body (a stick in two dimensions or a disk in three dimensions) is studied and it is proven that 
the system possesses a stationary velocity Voo which is also the limiting velocity and that the 
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asymptotic time behavior in approaching this Voo is power-law. More precisely, assuming the 
initial velocity Vq such that Wo — W) is positive and small, it was proven that: 

|14,-V(t)| a hg, (1.1) 

where d — 1, 2, 3 is the dimension of the physical space and C is a constant, depending on the 
medium and on the shape of the obstacle. This result, surprising for not being exponential, is 
due to re-collisions that can occur between gas particles and the body while it is accelerating. 

Similar model (HI ) have been studied where a stochastic kind of interaction between the 
gas and the body is assumed: when a particle of the medium hits the body it is absorbed and 
immediately stochastically re-emitted with a Maxwellian distribution centered around the body 
velocity, in this case the behaviour was found to be O () 

In f 1 H and |2j these models have been numerically studied, in particular they computed the 
dynamic with stochastic interaction carried out in JlJ for a disk subjected to an harmonic force 
confirming the analytical results (for analytical studies see 0 and references quoted therein), 
here they removed the hypothesis of initial velocity close to Vco, showing that this doesn't affect 
the dynamic. 

The question of whether the trend of the solution was connected with the simple shape 
chosen for the object in the above mentioned articles was firstly faced in 0, in the domain 
of elastic collisions, where it was studied the evolution of a general convex body, and it was 
confirmed the same power decay expressed in 0; this was mainly due to an important feature 
that the convex body shares with the first case of a simple disk, namely in both cases colliding 
gas particles bounce away from the body. 

In HOl the hypothesis of convexity was then removed and it was studied the case of a body 
with lateral barriers of finite length parallel to the axis of motion (a box-like object in two 
dimensions), it was shown that this simple kind of concavity significantly altered the asymptotic 
of the system; it was the first case in which there was a change in the power of decay: 0(f~ 3 ) 
and no change from two to three dimensions. The case of more general shape of the body was 
left out for further difficulties, we believe indeed that there are actually non trivial issues to be 
faced for these more general cases. 

To show this and in the in the spirit of getting a more realistic perspective in the study of 
friction problems, in this paper we studied the friction term itself acting on a body fixed at 
constant velocities but with a more general kind of concavity, for example the friction exerted 
by the air on a vehicle moving at constant velocity, ( as a practical image one can think to the 
air on the windshield of a car moving at constant speed) We show atypical features, namely the 
presence of a time dependent friction for constant velocities. 

In all previous cases in fact the object (the disk, the convex body or the box) fixed at any 
constant velocity experienced a constant friction. The reason being that a gas particle that hits 
a disk or a convex body moving at a fixed velocity cannot undergo a second collision with it, 
thus canceling the effect of recollisions and possible time dependencies in the friction, as far as 
the box shaped object in 1(101 , although recolisions can occur even when it is fixed at constant 
velocity, they will necessarily occur with the lateral barriers of the body which are parallel to 
the axis of motion, therefore the momentum transferred along the motion direction is zero, thus 
canceling again recollisions contributions. 

This is not true anymore if we consider the case of the present paper were a slightly more 
general kind of concavity is considered, namely a body with tilted lateral barriers (see Section 
[2] ), these trap particles which, bouncing inside, correlate through times different areas of the 
inner side of the body. The time dependence in the friction term is vanishing and we give an 
exact estimate of it. 

We then employ this result to prove that in the actual dynamic of the body (that is a body which 
is not fixed in velocity) there can't be a stationary velocity. The gas is assumed to be made of 
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free particles (see 0 on Knudsen gas) elastically interacting with the body and it is studied in 
the mean field approximation. 


2. The model 

We consider for sake of simplicity the two dimensional case. In particular we consider a 
symmetric angle-shaped body moving towards the x direction in a Vlasov gas, the body is 
constrained to stay with its base orthogonal to the x axis, with the center moving along the 
same axis and with the hollow base facing forwards, the thickness of the body is assumed 
negligible for sake of simplicity. The system is immersed in a perfect gas in equilibrium at 
temperature T and with constant density p, assumed in the mean field approximation, (that is 
the limit in which the mass of the particles goes to zero,while the number of particles per unit 
volume diverges, so that the mass density stays finite.) We will study the force exerted by the 
gas on to the body moving at a fixed velocity. We refer to our body as C(f) = C+(f) (J C_(t), 

C±(t) — {{x,y) E 1R 2 : x — X(f) + y cos(O); y — ±rj sin(0), rj E [0, L\} (2.1) 

where 6 E If, f) represents the angle between the x axis and the upper side of the body whose 
vertex is O(t) — (X(f),0). 

We refer to N(x), x E C(f) as the right unit normal vector to C(f), namely N(x) ■ x > 0; 
besides for further convenience we set fi = N (x) : x E C + (f), and ft — N(x) : x E C~ (f), in 
particular h — (sin(0), — cos(0)) and p — (sin(0),cos(0)) (see Figure 1). 

Let then f(x,v,t), (x, v) E IR 2 x ]R 2 be the mass density in the phase space of the gas 
particles. It evolves according to the free Vlasov equation: 


(df + v ■ \ , x )f{x,v,t) = 0, x<£C(t) (2.2) 

Together with eq. we consider the boundary conditions. They express conservation 
of density along trajectories with elastic reflection on C(f). Let then v be the velocity of a gas 
particle that hits the body at time t at the collision point R G C(f) and V(t) = X(f) the velocity 
of the body, denoting by Vn — v ■ N and by V\; = V(t)N ■ x, imposing elastic reflection in R, 
we have for the outgoing velocity w (for a derivation see the appendix): 

iv N = 2V N (t) - v N , w N± = v N± (2.3) 

where Vn ± — v — i> n N 

The boundary condition then reads as follows: 

f + {x,w,t) — f-(x,v,t) xeC(t) (2.4) 

where 

f±(x,v,t) — lim fix ± ev, v, t ± e). (2.5) 

e—M3+ 

Finally we give the initial state of the gas, assumed in thermal equilibrium through the 
Maxwell- Boltzmann distribution 

f{x,v,0) =fo{v 2 ) = p(^)e"^ 2 (2.6) 

7T 


with P=knT- 

The above equation for the gas are clearly coupled with those of the body immersed in it. 
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Figure 1: The body moving in the gas; the dashed line represents the trajectory of one gas particle hitting the body 


which are: 


fm = m, 

X(0) - 0, V(0) - % 

where E is the constant external force acting along the x-axis and 

L 

F(t) =4sin(0) J dp J dvjv 1 ■ n) 2 f-(X(t) + p cos(6),p sin(9),v,t) + 

0 l v' n<0 (2 8) 

—4sin(0) J dp j dv(v'■ n) 2 f-(X(t)+ pcos(6),psm(9),v,t) 

0 v'fi> 0 

is the action of the gas on the disk, being v' — v — V (t)x the velocity in the reference system of 
the body, and where fi — Nr for x G C+(f), in particular n — (sin(0), — cos(0)). 

We give here a derivation of eq. and | |2.8} . 

The contributions to F(f) coming from C+(f) and C_(f) are the same for the symmetry on the 
y axis of the gas distribution and the shape of the body, we thus consider the action of the gas 
on C+(f). Our body, while moving, is subjected to multiples collision with gas particles, if we 
write its total variation of momentum in an interval (f, t + At) along the x axis as A V(t) and 
the variation of momentum due to collisions with gas particles as AV co u(t), we have that (M is 
the mass of the body) 

MAV (t) — EAt + MAV co u{t). (2.9) 

In the following for any vector a we will use the notation: a„ — a ■ n ; — a ■ h j_ and for 

sake of simplicity we will denote V(t) simply as V. 


(2.7a) 

(2.7b) 

(2.7c) 
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After one collision at time t between a particle of position and speed (x,v) and C+(f) the 
change in momentum along x-axis is (see Appendix) (2 m/M){ v n — V n )h ■ x . 

The term AV co u(t) takes into account all the collisions happening during At, thus: 

A Vcoii (f) = ^ O v\ - V n )n ■ x + a (2.10) 

k 


where k labels all particles around the body that are hitting C+(f) within At, and a denotes 
terms o(At) 

Let Ax l Avl be a volume of the phase space of measure \Ax l AvI = AW, centered at the 
point ( x l , vl) and AN(x', vf t) the number of particles contained in it at time t, so that : 


^ E( v n - V " )n ■ A = TT E( v n ~ V » ' &A N(x l , vf t) 


M 


= M n*- v ») 


.. .. AN(x l , v ft) 

n- xm -—-AW, 

AW 


where i,j ranges over positions and velocities that will give rise to collision in (f, t + At). At 
this point letting AW —> 0 the mean field approximation,whose meaning was mentioned in the 
introduction, guarantees the convergence to a finite mass density, i.e : 


lim m 

AW—>0 


AN(xfvft) 

AW 


f(x, v,t), 


( 2 . 11 ) 


so that we arrive to: 


AV co u(t) — J ( v n — V n )h ■ xf(x, v, t)dxdv + a (2.12) 

n(Af) 

where Ci(At) is the (x,v) region of particles hitting C + (f) in ( t, t + At); for further convenience 
we split this integral into backward contribution to recoliision Cl 1 (At) and frontal contribution 
Ct R (A t): 


J — ( v n — V n )h ■ xf(x, v, t)dxdv + J —(v n — V n )fi-xf(x,v,t)dxdv. (2.13) 
n L (A t) n R (A t) 

These region are clearer if written through coordinates moving with the body, namely: 

i = x-0{t) (2.14) 

through this coordinates we can write C + (f) = {(£h,£_l) € 1R 2 : 0 < £j_ < L, = 0 } 

As far as ( f (At) is concerned, a collision happens only if v n — V n >0, besides the particle 
colliding are those contained in the rectangle: —{v n — V n )At < £ n < 0 up to order o(At). In this 
way we can compute the first integral as: 

L 0 

^ J dv J d£ ± J d£ n (v n -V n )n-xf(£ + 0(t),v,t) = 

V n > Vn 0 —[v n — V n )A t 

L 

= J dv J d£_ L {v n -V n ) 2 n-xf(€ + 0(t),v,t) + a 

V„>V„ 0 
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In the same way, as far as Q R (Af) is concerned, a collision happens only if v n — V„ < 0, besides 
the particle colliding are those contained in the rectangle: 0 < < (V n — v n )At up to order 

o(Af), so that for the second integral it yields: 

L (V„-z;„)A t 

^ J dv J d£ ± J d£ n (v n - V n )n-xf(€ + 0(t),v r t) = 

v„<V n 0 0 

L 

— ~f^At J dv J d£ ± (v n - V n ) 2 n-xf(£ + 0(t),v,t) + a 

v n <V„ o 


In the limit At —» 0 the integrand region becomes {(£,„ £j_) : 0 < £j_ < L, g n — 0 }, namely 
the C+(f), which moving back to the previous coordinates reads simply as 0 {(x,y) : x = 
X(f) + rj cos(0); y = ;/ sin(0), ?/ G [0, L]}, so we finally arrive to: 


lim AVc f ( ^ 
Af-)-0 At 


+2 h-xjdrj J dv( v n — V n ) 2 /_ (X(t) + rj cos(9), rj sin(0),u, t) + 

0 V„>V n 
L 

-2fi-xjdp J dv(v n — V n ) 2 /_(X(t) + ^cos(0),j;sin(0),u, t) 

0 v„<V n 


as we stated earlier the contribution to F(f) from the whole body is twice that on C + (f), besides 
noticing that h ■ x — sin (d) and v n — V„ — v' ■ ft concludes the derivation of eq.|2.7| and (2.8 . 


Now, thanks to boundary conditions we can compute the friction in terms of the initial 
mass density. Indeed let x(s,t,x,v), v(s, t, x, v) be the position and velocity of a particle at 
time s < t, that at time t occupies position x and velocity v ; conservation of mass imply that 
the P.d.f. stays constant along particles trajectories and in particular 

f(x, V, t ) = fo(x(0, t , X, u), 17 ( 0 , t , X, u)) (2.15) 


so that the problem of finding the gas distribution reduces to that of tracking the particles 
trajectories. 

Given the evolution of the body X(f), V ft), there is a unique backward time evolution 
leading to the initial position and velocity, (see Figure 1). Such backward evolution is free 
motion up to possible collision-times in which the particle hits the body. On these times we 
keep track of the particle displacement through condition | |2.4) . We proceed in this way until 
we reach the desired x{0,t,x,v) , v(0, t, x, v). At the end using the initial state of the gas 
distribution, eq.( 2 . 61 , we obtain 


— J J dv(v' ■ h) 2 e P v o — j dp J dv(v' ■ n) 2 e P v o ( 2 . 16 ) 

0 v' fi <0 0 v'ft >0 

where v 0 — v(0, t,X(t) + 7/cos (6),)] sin(0),u) and k = p(^). 

Note that in order to compute F(t) we need to evaluate v 0 and hence to know all the previous 
history (X(s), V (s), s < f}, on the other hand, if the light particle goes back without undergoing 
any collision, then vq — v and the friction term is easily computed 

Note also how the limit case "0-> | ", which gives back the normal stick in two dimensions, 
brings back the usual friction term, it suffices to note that under this limit v' ■ h = v x — V(t). 

1 see 2.1 
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3. Estimate of F v ( t ) 


In this section we will study the main problem, namely the friction exerted on a body fixed at 
(any) constant velocity. 

We start computing F v (f), the friction term for the body C(f) moving at a constant velocity, 
V ( t ) = V, in this case there can't be re collisions on the back of the body, though there can be 
on the front of it[^j therefore eq. 12.16 , in the case of a constant velocity, gets 


F V (t) 

4ksin(0) 



(3.1) 


with v' — v — Vx, for further convenience we will write V — Vx. 

In the first integral there is a subset of the region v' ■ n < 0 such that vq = v, that is particle 
not coming from a previous recollision, and a subset such that vq f= v, that is particles coming 
from a previous recollision, for these ones we now compute vo- 

As h is the normal to the right side of C+(f), we write p as the normal to the right side of 
C_ (t) , p — (sin(0),cos(0)) and pj_ = (— cos(0),sin(0)) as its upward normal; besides for any 
vector w we call iVp — w ■ p and Wp ± — w ■ p±. 

Having reduced to the C+(f) in the integral, particles coming from a previous collision can 
only come from C_(f); besides the condition 9 G [ J, ff] ensure at most one recollision for 
constant velocities of the body. In the reference system of the body the computation reduces to 
an elastic collision between a fixed tilted wall and a particle with velocity Vq — vq — V before 
and v' — v — V after the collision. We thus arrive to 


VQp — 2Vp Vp (3.2a) 

Vop = V p± (3.2b) 

This allow us to compute 

vl = v 2 3 -W v {v-V)-p (3.3) 

we stress again that this is the initial velocity of particles that underwent previous collision with 
C_ (f), we now compute the exact region in phase space that leads to recollisions. 

Let R(p) — (X(t) + 1 ] cos (8), ;/ sin(0)) be a point of C+(t) and Q — ( X(t ) + Lcos(0), — Lsin(9 )) 
be the extreme lower point of C_ (f), then 

$(*]) = N((t] + L) sin(0), (L - ?/) cos(0)) (3.4) 


is the vector such that 


^(f]) ■ RQ{rj) — 0 (3.5a) 

-0 ( 77 ) ■ h > 0 (3.5b) 

and it is clear that the subset leading to recollision is within the regiorj^j {?; t 1 R 2 : v' ■ h < 0 , 
v' ■ rf > 0 }, where v' — v — V. 

Now, to define the exact recollision region we follow backwards in time a particle hitting R{rj) 
with velocity v at time t. This particle will hit C_(t) if it covers in the p direction, a lenght 
greater than OR ■ p — rj sin(a) where oc — n — 29, which is to say ( v' ■ p)t > rj sin(20). 

2 This is indeed the main difference with all previous cases for constant velocities 

3 From now on we will simply write ip 
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Finally the exact recollision region reads 

TZg(p, t) — {v £ R 2 : v' ■ ft < 0, v' ■ xjr > 0, v' ■ p > y sin(20)} (3.6) 

Before proceeding we note that we can rewrite F v (f) in a convenient way: 

F v (t)=F 0 (V)+g(V,t) (3.7) 


with 


Fo(V) 

4 k sin(0) 



g(v,t) 

Ak sin(0) 


L 

j dp J dv(v r ■ n) 2 (e~P v o—e~P v2 ). 

0 v'fi< 0 


(3.8) 

(3.9) 


An important and new feature of this concave case is that even with constant velocity there 
is a time dependent term in the friction . Notice in fact that vo(x,v,t) — v for v IZgfp, t) 
hence leading to a null contribution in g(V, t), therefore we arrive to 


g(V,t) = 4fcsin(0) j dp j dvfv'-n) 2 e +v '> 2 [e i ^ v r^ v — 1 ], (3.10) 

o n e ( v ,t) 

Where we used eq.(3.3 1 because the integral is over the exact recollision region. [^] In this 
region particles had a previous collision thus VQp < Vp and through eq. 13.2 ■ it follows that 
v' ■ p > 0 . 

Moreover it is clear that 7tg(p,t) C Vff{p) where 


n^{p) = {i)£K 2 :»'4<0,))'i>0}, (3.11) 

this can be looked at as the recollision region for t —S- oo; these relations imply 

o <g(V,t)<g°°(V) (3.12) 

L 

g°°(V) — 4fcsin(0) f dp f dv'(v'-n) [ e ¥^K'P) _ i] ( 3 . 13 ) 

0 Uf( v ) 


this means that the concavity add a positive contribution g(V,t) to the friction and it is limited 
in time. 

In order to study the time dependence we write g(V,t) = g 00 (V) — A g(V, t ) with : 


*g(V,t)=g a> (V)-g(V,t) 


L 

Aksin(O) I dp f dv'(v'■n) 2 e-M v ' +v f[eW v p( v '-ft 
0 ATZ(t) 


1 ] 

(3.14) 


ATZ(t) — Tt™(p)\lZg(p, f) = {« £ R 2 : i) 1 ■ n < 0, v' ■ $ > 0, 0 < v' ■ p < j sin(20)} (3.15) 


We now perform the following rotation in the velocity space in order to normalize the integral 
4 Moreover since v' = v — V it yields dv' = dv. 
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region: 


v' — R(cp)w 



(3.16) 

(3.17) 


where tp = f — 0. We note that scalar products keep the same structure under R, that is for 
every vector q it yields 

v' ■ q — R((p)w ■ R((p)q w — w ■ q W/ (3.18) 

where q w — R(—<p)q; note also that p — (cos(f ),s\n(f )'). Through this change of variables, for 
which dw — dv r , we arrive to: 

L 

Ag(V,T) — Aksin(8)e~^ v j dq j dw(w-n l0 ) 2 e~^ w+2w ' Vw ^[e^ v P^ Wl ' > — 1] (3.19) 

0 AKw(T) 

A1Z W (T) — {re E 1R 2 : w ■ n w < 0, w ■ ip w > 0, 0 < W\ < Trj} (3.20) 

where tf w (i/) — N W (L — q cos (2d), —rj sin(26)), fi w — (— cos(20), — sin(20)) and T — sin ^ 2e \ 
besides w ■ fi w < 0 can be rewritten as w-\ n lV] + win W2 < 0 and, keeping in mind that n Wl < 0 
this leads to | "f 1 |tC] < u> 2 , similarly for w ■ ff> w > 0, leading us to : 

L Trj bwi 

Ag(V, T) — Aksin(8)e~P v " J dij j dw\ J dw 2 f(w\,W 2 ) (3.21) 

0 0 cixv\ 

f(w 1 , W 2 ) - (w ■ n w ) 2 e-^ w2+2w V ^ [e 4 P v rM _ i] (3.22) 


where 

a — I ~~ I — tan (20 ~ y) 

W-IV2 ^ 

, _ | ipw! | _ L - ?/cos(2 9) 

ipy 2 q sin ( 20 ) 

we incidentally observe that b > tan( 0 ), which we will be of use later. 


(3.23) 

(3.24) 


We now give an exact estimate in time of A g(V,T). In the sequel C and C n will denote 
positive c onstan ts possibly depending on p, 8, V and L. 

From eq. (3.14 1, using that e~^ v ' +v 'i < 1 , (v’ ■ ft,) 2 < \v’\ 2 and passing to the variables "w" 
we arrive to 

L 

Ag(V,t) <C f drj f dw w 2 [e c llUl - 1], (3.25) 

0 A Tlw 


clearly aq < Trj and w 2 < \(Tq,bTq)\ 2 < CT 2 since, from 1 3.24 1 , bi) is boimded, moreover 
\ATl w \ < CT 2 , this leads us to 


and then 


A g(V,T) < CT 4 (e c > T - 1) 


lim 

T—>0 


Ag(VJ) 

p5 


< c. 


(3.26) 

(3.27) 
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For the lower bound we start again from eq. 
leading us to 

L 

A g(V,t)>C J dp J dvfv'■h) 2 e~^ v '\ 2 + 2 \ v 'W v \)[ e i l i M v '-P) - 1} (3.28) 

L/2 AK(t) 


3.14 1 and use that \v' + V \ 2 < (|t/| + | V\ 


Now passing to the variables "w" and using again that w 2 < CT 2 , 

we have that w 2 + 21 w\ V < CT 2 + C\T < CfT, the last inequality holding for T small enough; 
we thus arrive to 

L Tij b w 1 

Ag(V,t)>C J dp J dw\ I du? 2 (w ■ n w ) 2 e~ ClT [e C2Wl — 1 ] > (3.29) 

L/2 0 a wi 

L Tt] b w 1 

> C J dp J dw\ J dw 2 T 2 e~ ClT [e ClT — 1 ] > (3.30) 

L/2 Ty/2 «e 

> CT 4 e" ClT [e C2T - 1], (3.31) 

in the last inequalities we shrunk the integral region in: a £ W\ < zv 2 < bw where 

a e — tan (20 — tz/2 + e) can be chosen satisfying: tan (20 — n/2) — a < a £ < tan( 0 ) < b, with 

0 <e< 7 r /2 — 0;in this region it holds: 


(w n w ) — \n W2 \(azv 1 - zv 2 ) < - \n W2 \{a e - a)zv 1 < -CT 


leading to ( w ■ n w ) 2 > CT 2 which was used in 13.30 . Therefore we arrived to 


lim 

T—>0 


A g(V,T) 

p5 


> c. 


(3.32) 


(3.33) 


Now, changing back to t — sm ^\ together with the fact that A g(V, t) is bounded, leads us to 
the desired estimate for Ag(V,t), namely 


C? < Ag(V, t ) < Cm ,^ Vf > 0 


(1 + t ) 5 


(1 +t) 5 


and finally 


- (TT7)5 ^ pV W - ^ - (TTfp vt ~ °' 


(3.34) 


(3.35) 


4. Absence of a stationary solution 


We now show how having computed the fricion term in the previous problem also allow us 
to prove that the actual dynamic ( |2.7| ) of such a concave body (no more fixed at any constant 
velocity) has no stationary solution. For the previous cases of the disk, the convex body and 
the cylinder there is a stationary solution V(t) = V that can be found directly from eq. 1 2.7 1 
namely : 


0 = E-F V (t), Vt > 0 
V(0) - V 


(4.1) 

(4.2) 
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where F v (f) is the friction term computed for the constant velocity V; in those cases (see 0, 
PI and 111 01 ) it yields F v (t) = Fq (V ) Vf > 0 because recolision terms are absent for constant 
velocities, the stationary solution is therefore the one such that E — Fq(V) = 0 which was called 
Voo and was also the limiting velocity of the body 

In the present case of concavity the possible stationary velocity has to be such that 


0 = E-F v (t) =E-F 0 (V)-g{V,t), Vt > 0 
1/(0) = V. 


(4.3) 

(4.4) 


This implies the following necessary condition on the stationary velocity: 


dg(V, t) Q 


Vf > 0 


(4.5) 


and through the previous analysis we can compute this derivative from 1 3.21} in the following 
way (we use for our convenience the variable employed earlier T = sin (20) 


dg(V,T) 

dT 


L br/T 

T ^ = 4A: sin(6)e~P v2 j dp // j dw 2 f(pT,w 2 ) 

0 arjT 


(4.6) 


We immediately observe that (4.3 ■ is not satisfied by V = 0 since Fo(0) — g(0, f) = 0, while for 
1? > 0 we observe by immediate inspection that f( i]T, iv 2 ) >0 VV > 0, thus 


dg(V,T) 

dT 


< 0 


W > 0, T > 0 


(4.7) 


which exclude the possibility of a stationary solution. 


5. Comments 

In this work we studied what we believe to be the main feature of general kinds of concavities 
in their interaction with a fluid, namely tilted walls that correlate the body with itself via 
the bouncing of gas particles, we also showed how this affects even the simple friction for 
constant velocities, in particular the presence of a time dependence in it. As it was reasonable 
to expect, not only this contribution is bounded but also vanishing in time trough the estimate 
we computed. 

We studied the two dimensional case in order to focus on the main feature of the problem, 
avoiding long technical details; following the reasoning of the present work it is evident that the 
only tweak between two and three dimensions is the exact power at which the time dependence 
vanishes, the rest being the same. 

We also stress that these results, as those of the main articles to which we refer, strongly 
depend on the choice of a perfect gas as the mean exerting friction, recollisions in fact come 
essentially from the long memory effect of gas particles in the setting of a Vlasov gas. 

We end this section observing that the absence of a stationary solution doesn't exclude a limiting 
velocity in the dynamic, which is instead expected, in fact we showed that the time dependence 
of this new friction term is vanishing reasonably fast in time and through heuristic reasoning 
we expect the limit velocity to be a Wo satisfying 

F — Fo(Voo) — g°°(Voo) = 0 (5.1) 

which is a correction of the previous Voo ( see j4j, f5| and IflOl ) due to the non vanishing com¬ 
ponent of the concavity friction, thus depending also on the tilting parameter 0, in particular, as 


11 












Time dependent friction in a free gas 


it is clear from the foregoing analysis, Voo —> Voo for 6 —> 7t/2. 

The actual dynamic for this kind of general concavity is still to be studied, but this work 
together with [10] show how the shape of the body, in particular a meaningful concavity, is 
essential in the interaction with the medium it is moving in. 


A. Appendix 

We derive here collision conditions (2.3) and (2.4). In what follows we denote by M and V the 
mass and velocity of the body and by m and v mass and velocity of a particle which will be 
assumed to collide elastically with the body. 

We consider that the particle hits the body at x G C(f) at time t and define its pre(post)collisional 
velocity as v (?/), projecting conservation of momentum equation on N and N L direction, 
together with the conservation of kinetic energy imply for the body: 

2.771 

V N = V N + M _|_ m (vn ~ V N ) (A.l) 

V^ ± = V Nx (A.2) 

where Vf and V N are post-collisional velocities along the N direction. 

While for the particle: 


v n ~ Vn 


M — m 
M + m 


(v N ~ V N ) 


v 'n ± = vn ± 


(A3) 

(A.4) 


Being M » m, we have 


v' N ~ 2V n - v N 

V N ± = ^N ± - 

and for the body: 

2tti 

v n — v n + ^ (p N - v N ) 

Therefore 

2 th /v 

V —V — (vfi — V n )N = — (v N - V n )N 
hence the momentum change along the x direction is 

2/77 

— (t’N - V n )N -X. 


(A.5a) 
(A.5b) 


(A.6) 

(A .7) 


(A.8) 
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